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Abstract 



The fully developed turbulence with weak anisotropy is investigated by 
means of renormalization group approach (RG) and double expansion regu- 
larization for dimensions d > 2. Some modification of the standard minimal 
substraction scheme has been used to analyze stability of the Kolmogorov scal- 
ing regime which is governed by the renormalization group fixed point. This 
fixed point is unstable at d = 2; thus, the infinitesimally weak anisotropy de- 
stroyes above scaling regime in two-dimensional space. The restoration of the 
stability of this fixed point, under transition from d = 2 to d = 3, has been 
' demonstrated at borderline dimension 2 < d c < 3. The results are in qual- 

itative agreement with ones obtained recently in the framework of the usual 
analytical regularization scheme. 

1 Introduction 

A traditional approach to the description of fully developed turbulence is based on 
the stochastic Navier-Stokes equation Jl|. The complexity of this equation leads to 
great difficulties which do not allow one to solve it even in the simplest case when 
one assumes the isotropy of the system under consideration. On the other hand, the 
isotropic turbulence is almost delusion and if exists is still rather rare. Therefore, 
if one wants to model more or less realistic developed turbulence, one is pushed to 
consider anisotropically forced turbulence rather than isotropic one. This, of course, 
rapidly increases complexity of the corresponding differential equation which itself 
has to involve the part responsible for description of the anisotropy. An exact solution 
of the stochastic Navier-Stokes equation does not exist and one is forced to find out 
some convenient methods to touch the problem at least step by step. 

A suitable and also powerful tool in the theory of developed turbulence is the 
so-called renormalization group (RG) methodQ Over the last two decades the RG 

1 Here we consider the quantum-field renormalization group approach Q rather than the Wilson 
renormalization group technique H 
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technique has widely been used in this field of science and gives answers to some 
principal questions (e.g. the fundamental description of the infra-red (IR) scale 
invariance) and is also useful for calculations of many universal parameters (e.g. 
critical dimensions of the fields and their gradiens etc.). A detailed survey of this 
questions can be found in Refs. [|, [5| and Refs. therein. 

In early papers, the RG approach has been applied only to the isotropic mod- 
els of developed turbulence. However, the method can also be used (with some 
modifications) in the theory of anisotropic developed turbulence. A crucial question 
immediately arises here, whether the principal properties of the isotropic case and 
anisotropic one are the same at least at the qualitative level. If they are, then it 
is possible to consider the isotropic first step in the investigation of real 

systems. On this way of transition from the isotropic developed turbulence into the 
anisotropic one we have to learn whether the scaling regime does remain stable under 
this transition. That means, whether the stable fixed points of the RG equations 
remain stable under the influence of anisotropy. 

Over the last decade a few papers have appeared in which the above question has 
been considered. In some cases it has been found out that stability really takes place 
(see,e.g. |5], [7|]). On the other hand, existence of systems without such a stability has 
been proved too. As has been shown in Ref. ||, in the anisotropic magnetohydro- 
dynamic developed turbulence a stable regime generally does not exist. In |7], [J, 
d-dimensional models with d > 2 were investigated for two cases: weak anisotropy 
@ and strong one ||, and it has been shown that the stability of the isotropic fixed 
point is lost for dimensions d < d c = 2.68. It has also been shown that stability 
of the fixed point even for dimension d — 3 takes place only for sufficiently weak 
anisotropy. The only problem in these investigations is that it is impossible to use 
them in the case d = 2 because new ultra-violet (UV) divergences appear in the 
Green functions when one considers d = 2, and they were not taken into account in 

In [ 3D | , a correct treatment of the two-dimensional isotropic turbulence has been 
given. The correctness in the renormalization procedure has been reached by in- 
troduction into the model a new local term (with a new coupling constant) which 
allows one to remove additional UV-divergences. From this point of view, the results 
obtained earlier for anisotropical developed turbulence presented in [[Ll]] and based 
on the paper JT2] (the results of the last paper are in conflict with |r0|| ) cannot be 
considered as correct because they are inconsistent with the basic requirement of the 
UV-renormalization, namely with the requirement of the localness of the countert- 
erms 



mm. 



The authors of the recent paper \TE] have used the double-expansion procedure 
introduced in |nj (this procedure is a combination of the well-known Wilson dimen- 
sional regularization procedure and the analytical one) and the minimal substraction 
(MS) scheme [[TjJ for investigation of developed turbulence with weak anisotropy for 
d = 2. In such a perturbative approach the deviation of the spatial dimension from 
d = 2, S = {d — 2)/2, and the deviation of the exponent of the powerlike correlation 



2 



function of random forcing from their critical values, e, play the role of the expansion 
parameters. 

The main result of the paper was the conclusion that the two-dimensional fixed 
point is not stable under weak anisotropy. It means that 2d turbulence is very 
sensitive to the anisotropy and no stable scaling regimes exist in this case. In the case 
d = 3, for both the isotropic turbulence and anisotropic one, as it has been mentioned 
above, existence of the stable fixed point, which governs the Kolmogorov asymptotic 
regime, has been established by means of the RG approach by using the analytical 
regularization procedure |], 0, |9j . One can make analytical continuation from d = 2 
to the three-dimensional turbulence (in the same sense as in the theory of critical 
phenomena) and verify whether the stability of the fixed point (or, equivalently, 
stability of the Kolmogorov scaling regime) is restored. From the analysis made in 
Ref. Jy| it follows that it is impossible to restore the stable regime by transition from 
dimension d — 2 to d — 3. We suppose that the main reason for the above described 
discrepancy is related to the straightforward application of the standard MS scheme. 
In the standard MS scheme one works with the purely divergent part of the Green 
functions only, and in concrete calculations its dependence on the space dimension 
d resulting from the tensor nature of these Green functions is neglected (see Sect. 3). 
In the case of isotropic models, the stability of the fixed points is independent of 
dimension d. However, in anisotropic models the stability of fixed points depends on 
the dimension d, and consideration of the tensor structure of Feynman graphs in the 
analysis of their divergences becomes important. 

In present paper we suggest applying modified MS scheme in which we keep the 
d— dependence of UV-divergences of graphs. We affirm that after such a modification 
the (i-dependence is correctly taken into account and can be used in investigation 
of whether it is possible to restore the stability of the anisotropic developed turbu- 
lence for some dimension d c when going from a two-dimensional system to a three- 
dimensional one. In the limit of infinitesimally weak anisotropy for the physically 
most reasonable value of e = 2, the value of the borderline dimension is d c = 2.44. 
Below the borderline dimension, the stable regime of the fixed point of the isotropic 
developed turbulence is lost by influence of weak anisotropy. 

It has to be mentioned that a similar idea of "geometric factor" was used in 
Ref. ]17|] in RG analysis of the Burgers-Kardar-Parisi-Zhang equation but the reason 
to keep the d— dependence of divergent parts of the graphs was to take correctly into 
account the finite part of one-loop Feynman diagrams in the two-loop approximation. 
In the present paper, we shall not discuss it in detail because the critical analysis of 
the results obtained in [17] was given in 



The paper is organized as follows. In Section 2 we give the quantum field 
functional formulation of the problem of the fully developed turbulence with weak 
anisotropy. The RG analysis is given in Section 3 when we discuss the stability of 
the fixed point obtained under weak anisotropy. In Section 4 we discuss our results. 
Appendix I contains expressions for the divergent parts of the important graphs. 
At the end, Appendix II contains analytical expressions for the fixed point and the 
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equation which describes its stability in the limit of the weak anisotropy. 

2 Description of the Model. UV-divergencies 

In this section we give the description of the model. As has already been discussed 
in the previous section, we work with fully developed turbulence and assume weak 
anisotropy of the system. It means that the parameters that describe deviations from 
the fully isotropic case are sufficiently small and allow one to forget about corrections 
of higher degrees (than linear) which are made by them. 

In the statistical theory of anisotropically developed turbulence the turbulent 
flow can be described by a random velocity field v(x, t) and its evolution is given by 
the randomly forced Navier-Stokes equation 

*L + @.$)0- V)A v-fA = f i (i) 

where we assume incompressibility of the fluid, which is mathematically given by the 
well-known conditions V • v — and V ■ / = 0. In eq.(|l|) the parameter Uq is the 
kinematic viscosity (hereafter all parameters with subscript denote bare parameters 
of unrenormalized theory, see below), the term f A is related to anisotropy and will 
be specified later. The large-scale random force per unit mass / is assumed to have 
Gaussian statistics defined by the averages 

(fi) = 0, (M£i,t)fj(£ 2 ,t)) = D l3 [x x - x 2 ,h - t 2 ). (2) 

The two-point correlation matrix 

r d d k ~ 

Dij{x, t) = 5(t) J ——D^k) exp(ik ■ x) (3) 

is convenient to be parameterized in the following way || ||: 

D l0 {k) = g Q vlk i - d - 2 %l + a 10 e k )P t ,(k) + a 20 Ri,(k)] , (4) 

where a vector k is the wave vector, d is the dimension of the space (in our case: 
2 < d), e > is a dimensionless parameter of the model. If the dimension of the 
system is taken d > 2, then the physical value of this parameter is e = 2 (the so- 
called energy pumping regime). The situation is more complicated when d = 2. In 
this case the new integrals of motion arise, namely the enstrophy, and all its powers 
(for details see Ref. ||19|| ) which leads to ambiguity in determination of the inertial 
range and this freedom is given in the RG method by the value of the parameter 
e. The value e = 3 corresponds to the so-called enstrophy pumping regime. This 
problem of uncertainty cannot be solved in the framework of the RG technique. On 
the other hand, the value of e is not important for stability of the fixed point when 
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d = 2. Thus, it is not important, from our point of view, what is the value of e in 
the case d = 2. Its value e = corresponds to a logarithmic perturbation theory 
for calculation of Green functions when g , which plays the role of a bare coupling 
constant of the model, becomes dimensionless. The problem of continuation from 



e = to physical values has been discussed in |2T) |. The (d x <i)-matrices and Rij 
are the transverse projection operators and in the wave-number space are defined by 
the relations 

Pij(k) = % - ' Ri i$) = (ni - ikj^j (nj -ZkjJ , (5) 

where is given by the equation = k ■ n/k . In eq.(|5|) the unit vector n specifies 
the direction of the anisotropy axis. The tensor given by eq.(|3]) is the most 
general form with respect to the condition of incompressibility of the system under 
consideration and contains two dimensionless free parameters ctio and a 2 Q. From 
the positiveness of the correlator tensor one immediately gets restrictions on 
the above parameters, namely «io > — 1 and a 2 Q > 0. In what follows we assume 
that these parameters are small enough and generate only small deviations from the 
isotropy case. 

Using the well-known Martin-Siggia-Rose formalism of the stochastic quantiza- 
tion [^T| 22, |24| one can transform the stochastic problem (P with the correlator 



(|^) into the quantum field model of the fields v and v\ Here V' is an independent of 
the v auxiliary incompressible field which we have to introduce when transforming 
the stochastic problem into the functional form. 

The action of the fields v and w is given in the form 

S = - J d d xidt 1 d d x 2 dt 2 vl(xi,ti)Dij(xi - x 2 -,ti -t 2 )v , j (x2-,t 2 ) 

+ J d d xdt [v< (x, t) [-d t v - {v ■ V)v + u V 2 v + f A ] {x, t) } . (6) 

The functional formulation gives the possibility of using the quantum field theory 
methods including the RG technique to solve the problem. By means of the RG 
approach it is possible to extract large-scale asymptotic behaviour of the correlation 
functions after an appropriate renormalization procedure which is needed to remove 
U V- divergences . 

Now we can return back to give an explicit form of the anisotropic dissipative term 
f A . When d > 2 the UV-diver gences are only present in the one-particle-irreducible 
Green function < v"u >. To remove them, one needs to introduce into the action 
in addition to the counterterm v'V 2 v (the only counterterm needed in the isotropic 
model) the following ones v'{n- V) 2 w, (n ■ v^V 2 (n ■ v) and (n-v')(n- V) 2 {n-v). These 
additional terms are needed to remove divergences related to anisotropic structures. 
In this case (d > 2), one can use the above action (0) with ([|) to solve the anisotropic 
turbulent problem. Therefore, in order to arrive at the multiplicative renormalizable 
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model, we have to take the term f A in the form 

f A = "o [xio(n • V) 2 v + X2o^V 2 (n ■ v) + Xmn(n ■ V) 2 (n ■ v)\ . (7) 

Bare parameters Xio> X20 and X30 characterize the weight of the individual structures 
in®. 

A more complicated situation arises in the specific case d = 2 where new diver- 
gences appear. They are related to the 1-irreducible Green function < v"tJ' > which 
is finite when d > 2. Here one comes to a problem how to remove these divergences 
because the term in our action, which contains a structure of this type is nonlo- 
cal, namely v'k 4 ~ d ~ 2e v'. The only correct way of solving the above problem is to 
introduce into the action a new local term of the form t?'V 2, y' (isotropic case) [TO. 
In the anisotropic case, we have to introduce additional counterterms v'(n ■ V) 2 v', 
(n ■ v')V 2 {n ■ v>) and (ft ■ v>)(n ■ V) 2 (n • v>). In [|TI], [T2[ a double-expansion method 
with a simultaneous deviation 25 = d — 2 from the spatial dimension d = 2 and 
also a deviation e from the k 2 form of the forcing pair correlation function propor- 
tional to k 2 ~ 2S ~ 2e was proposed. We shall follow the formulation founded on the 
two-expansion parameters in the present paper. 

In this case, the kernel (f|) corresponding to the correlation matrix D^ixi — 
x 2 ,t 2 — t\) in action (H is replaced by the expression 

+ 92^ k \ l + ^l)^{k) + (a m + a^QR ij {k)}. (8) 

Here and are given by relations (^), g 2 o, 0:30, «4o an d 050 are new parameters 
of the model, and the parameter go in eq. (f§) is now renamed as gio- One can see that 
in such a formulation the counterterm v , V 2 v' and all anisotropic terms can be taken 
into account by renormalization of the coupling constant g 2 o, and the parameters 
a 30 , a 4 o and a 50 . 

The action (§) with the kernel Dij(k) (§) is given in the form convenient for 
realization of the quantum field perturbation analysis with the standard Feynman 
diagram technique. From the quadratic part of the action one obtains the matrix of 
bare propagators (in the wave-number - frequency representation) 



— =< ViVj > = AjJ(k, u k ), 

+- =< vtv'j > = A^J' (k, w fc ), 

+- = < viv' j > Q =Av?'(k,u k ) = 0, 



where 



K: 
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+ 



1 



K x {K 2 + K{l-tl)) 



KK 3 + K(K 3 + K 4 (l-e k )) 



( Kl +K(i-e k )) 



K 4 



Rij 



A vv 



(k,LU k ) 



K 



K 2 (K 2 + K(l-e k )) 



-Ra 



ij ■ 



(9) 



with 



K 



iu k + u k 2 + u xw(n ■ k) 2 , 
-iu k + v k 2 + voXw(n ■ k) 2 
-g 10 ^k 2 - 2S - 2e (l + a l0 e k )- 

,3;„2-2<5-2<e. 



g 20 is$k 2 (l + awtl) 



"oX2ok 2 + v Q Xw{n • k) 2 . 



(10) 



The propagators are written in the form suitable also for strong anisotropy when 
the parameters are not small. In the case of weak anisotropy, it is possible to 
make the expansion and work only with linear terms with respect to all parameters 
which characterize anisotropy. The interaction vertex in our model is given by the 
expression 

j 



Viji = i(kj5u + kidi 



Here, the wave vector k corresponds to the field v\ Now one can use the above 
introduced Feynman rules for computation of all needed graphs. 

3 RG-analysis and Stability of the Fixed Point 

Using the standard analysis of quantum field theory (see e.g. [|], ||, [13], [14|]), one 



can find out that the UV divergences of one-particle-irreducible Green functions 
< vv> >jji and < vv >m are quadratic in the wave vector. The last one takes 
place only in the case when dimension of the space is two. All terms needed for 
removing the divergences are included in the action (0) with (0) and kernel (|8|). 
This leads to the fact that our model is multiplicatively renormalizable. Thus, one 
can immediately write down the renormalized action in wave-number - frequency 
representation with V — > ik, d t — > —ioj k (all needed integrations and summations are 
assumed) 

1 



S R (v,v) = -v\ 



gi u 3 ^k 2 - 2S - 2 " ((l + otig) P tJ + a 2 R i3 
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+ g 2 v 3 fi- 25 k 2 ((Z 5 + Z 6 a 3 Zl) Pij + (z 7 a 4 + Z 8 a 5 £ 2 ) Rij 



+ vi 



(iu k - Z x vk 2 )Pij - vk 2 {z 2 X\Ck P ij + (^3X2 + Z A X3^l) R%. 



+ \ v l v j v i V iji , (11) 

where \x is a scale setting parameter with the same canonical dimension as the wave 
number. Quantities g^ Xi, 013, «4, «5 and v are the renormalized counterparts of bare 
ones and Z { are renormalization constants which are expressed via the UV divergent 
parts of the functions < vv >jr and < vv >ir. Their general form in one loop 
approximation is 

Zi = 1 - Fi Poles? 6 . (12) 

In the standard MS scheme the amplitudes Fj are only some functions of 
Xi, o?3, «4, «5 and are independent of d and e. The terms Poles*' 6 are given by 
linear combinations of the poles \, | and (for 5 — > 0, e — > 0). The amplitudes 
Fj = F/F 2 are a product of two multipliers F/,F 2 . One of them, say, F/ is a 
multiplier originating from the divergent part of the Feynman diagrams, and the 
second one Ff is connected only with the tensor nature of the diagrams. We explain 
that using the following simple example. Consider an UV-divergent integral 

T(lr ^ - r, „ h h [j* n 1 mm™ Sijqiqm + Sgqjqm + Sjiqiqm 

i{K,n) = riiUjKiKm j a q + m2 y +2 &\ —± 3^ > 

(summations over repeated indices are implied) where m is an infrared mass. It can 
be simplified in the following way: 

J(k,n) = nirijkikraSijira dq - m2 y +25 , 



where 



q 25 T{6 + 1)T(6) 



dq 



2 



2{q 2 + m 2 y+ 25 2m 25 T{25 + 1)' 

and Sd = 2ix d l 2 /T(d/2) is the surface of unit the d-dimensional sphere. The purely 
UV divergent part manifests itself as the pole in 25 = d — 2; therefore, we find 

UV div. part of / = -^(F 2 A; 2 + F 2 2 (nk) 2 ), 

2o 

where F 2 = F 2 /2 = (1 - d)S d /3d(d + 2). 
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In the standard MS scheme one puts d = 2 in F 2 , F 2 ; therefore the d-dependence 
of these multipliers is ignored. For the theories with vector fields and, consequently, 
with tensor diagrams, where the sign of values of fixed points and/or their stabil- 
ity depend on the dimension d, the procedure, which eliminates the dependence of 
multipliers of the type F 2 , F 2 2 on d, is not completely correct because one is not able 
to control the stability of the fixed point when drives to d — 3. In the analysis of 
Feynman diagrams we propose to slightly modify the MS scheme in such a way that 
we keep the d-dependence of F in fll2|) . The following calculations of RG functions 
(/3— functions and anomalous dimensions) allow one to arrive at the results which 
are in qualitative agreement with the results obtained recently in the framework of 
the simple analytical regularization scheme ||, i.e. we obtain the fixed point which 
is not stable for d = 2, but whose stability is restored for a borderline dimension 
2 < d c < 3. 

The transition from the action (§) to the renormalized one (pTTJ ) is given by the 
introduction of the following renormalization constants Z: 



u = vZ u , 5-io = gifi 2e Z 9l , 5-20 = 92[i M Z 92 , Xio = XiZ 



-26, 



Xi ' 



a 



(i+2)0 



Oti+2Z ai+2 , (13) 



where i = 1,2, 3. By comparison of the corresponding terms in the action (p]) 
with definitions of the renormalization constants Z for the parameters (|T3"D, one can 
immediately write down relations between them. Namely, we have 



Z, 

z 



Hi 



J 9'2 



V-3 

Zc,Z 1 3 , 
Z\+iZ l 1 , 



(14) 



where again i — 1, 2, 3. 

In the one-loop approximation, divergent one-irreducible Green functions 
< w >m and < vv >m are represented by the Feynman graphs 



< vv > IR ~- 




< vv > IR = 




(15) 



The divergent parts of these diagrams T v ""' , T v ' v have the structure 



X 



+ 



+ 



9t 



Ae + 25 
2e 

A_ 

-25 



(ai5ijk 2 + a 2 5y(n ■ k) 2 + a 3 riinjk 2 + a$iinj{n ■ k)' 



bidijk + b2Sij(n ■ k) + b 3 ni7ijk + ^rijn^n • k)' 
{c\5ijk 2 + C25ij(n ■ k) 2 + c^nirijk 2 + c^riin^n ■ k)' 
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X 



9i 
2e 



(didijk 2 + d 2 5ij(n ■ k) 2 + d^niUjk 2 + d^riiUj (n ■ k) 2 ^j 
+ — j-z (eidijk 2 + e 2 5ij(n • k) 2 + e^riirijk 2 + e^niUj{n ■ k) 2 

£0 



(16) 



where parameter A and functions Oj, foj, Cj, di and are given in the Appendix I 
(i = 1...4). The counterterms are built up from these divergent parts which lead to 
the following equations for renormalization constants: 



Zi 



z x . 



1 A (9i, 92 
1 - U dl - 25 61 
, A (91, 92 
Xi V2e 26 



9i, 



92 



2 U 2 4e + 25 2e 25 , 



1 



.4 



9\ o-i+i . 9i, 

+ 7T°i+l 



1oi i+2 \g 2 ±e + 26 2e 
1,2,3. 



92 

26 



Ci+l 



(17) 



From these expressions one can define the corresponding anomalous dimensions 
7i = /ic^ In Zi for all renormalization constants Z\ (logarithmic derivative /i<9 M is 
taken at fixed values of all bare parameters). The /3-functions for all invariant charges 
(running coupling constants gi, g 2 , and parameters x%i a i+2) where i = 1,2,3 are 
given by the following relations: /3 gi = fid^i (i = 1, 2), (3 Xi = fid^Xi and (3 ai+2 = 
fid^ai + 2 (i = 1,2,3). Now using equations ( |14|) and definitions given above one can 
immediately write the /3-functions in the form 



-gi(2e + j gi ) = # 1 (-2e + 3 Tl ), 

^(25 - 7 92 ) = #2(25 + 37! - 75) 
-Xilxi = -Xiili+i - 7i) , 

-a i+ 2la l+2 = -aj+ 2 (7i+5 - 75) ■ 



1,2,3, 



where 



7i 
Ji+i 

75 



A{gidi + g 2 e x ) , 

— {sxdi + x + g2e i+ i) , 
Xi 

A (9l 

— — —a-i + gih + g 2 ci 

2 \92 



A 



'9l 



2a i+2 \g 2 



a i+l + 9lbi+l + 92<k+l 



1,2,3 



(19) 
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Figure 1: Dependence of the borderline dimension d c on the parameter e and for 
concrete values of a% and a 2 - 

By substitution of the anomalous dimensions 7« ( |I13D into the expressions for the 
/5-functions one obtains 



An 



gi(-2e + 3A(pidi + # 2 ei)) , 

92 



A ( g 2 

25 + 3A(g 1 d 1 + g- 2 ei) + 77 — a x + + g 2 c 1 

2 V^2 



A 



Jl a i+ i + g\b i+1 + g 2 c i+ i) + a i+2 f —a x + gih + g 2 c x 



i = 1,2,3. (20) 
The fixed point of the RG-equations is defined by the system of eight equations 

Pc(C.)=0, (21) 
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where we denote C = {gi, g 2 , x%, a i+2}, i = 1,2,3 and C* is the corresponding value 
for the fixed point. The IR stability of the fixed point is determined by the positive 
real parts of the eigenvalues of the matrix 

I, m=l, ...8. (22) 

c=c* 

Now we have all necessary tools at hand to investigate the fixed points and their 
stability. In the isotropic case all parameters which are connected with the anisotropy 
are equal to zero, and one can immediately find the Kolmogorov fixed point, namely: 

1_ 8(2 + d)e(2e - 3d(5 + e) + d 2 (35 + 2e)) 
A 9(-l + dfd{l + d){5 + e) ' 

1 8(-4 - 2d + 2d 2 + d 3 )e 2 

A 9(-l + d) 3 d(l + d)(S + e) ' ( j 




^2* = 



12 



where 5 = (d— 2)/2 and the corresponding matrix has the following eigenvalues: 

1 

6d(d — I) 



Ai j2 = — — -\m{d- 1) + 4e(2-3d + 2d 2 



± 



(6d<5(l - d) - 4e(2 - 3d + 2d 2 )) 2 



- 12d(d - l)e(12d<5(d - 1) + 4e(2 - 3d + 2d 2 )) 



(24) 



By a detail analysis of these eigenvalues we know that in the interesting region of 
parameters, namely e > and 5 > (it corresponds to d > 2) the above computed 
fixed point is stable. In the limit d = 2, this fixed point is in agreement with that 



given in |T(| |15 | 



When one considers the weak anisotropy case the situation becomes more com- 
plicated because of necessity to use all system of /5-functions if one wants to analyze 
the stability of the fixed point. It is also possible to find analytical expressions for 
the fixed point in this more complicated case because in the weak anisotropy limit 
it is enough to calculate linear corrections of a,\ and a 2 to all the quantities (see 
Appendix II). 

To investigate the stability of the fixed point it is necessary to apply it in the 
matrix (p2|) . Analysis of this matrix shows us that it can be written in the block- 
diagonal form: (6 x 6) (2 x 2). The (2 x 2) part is given by the /3-functions of the 
parameters and X3 an d, namely, this block is responsible for the existence of 
the borderline dimension d c because one of its eigenvalues, say Ai(e, d, a±, ct 2 ), has 
a solution d c e (2,3) of the equation Ai(e, d c , a±, a> 2 ) = for the defined values of 
e, ai,a2- The following procedure has been used to find the fixed point: First we 
have used the isotropic solution to g\ and g 2 to compute the expressions for c*i +2 and 
Xi, i = 1,2,3. From equations /3 a& = and f3 X3 = one can immediately find that 
«5* = and X3* — 0. After this we can calculate expressions for the fixed point of 
the parameters a i+2 and ^,1 = 1,2. At the end, we come back to the equations for 
gi and g 2 , namely /3 gi = and (3 92 = 0, and find linear corrections of a,\ and ot 2 to the 
fixed point. The corresponding expressions for the fixed point and the corresponding 
eigenvalue of the stability matrix responsible for instability are given in Appendix 
II. 

From a numerical analysis of the stability matrix one can find that in some region 
of space dimensions d the stability is lost by the influence of the weak anisotropy. 
On the other hand, the borderline dimension d c arises when going from dimension 
d = 2 to d = 3. For the energy pumping regime (e = 2) we found the critical 
dimension d c = 2.44. This value corresponds to ct\ = a 2 = 0. This is the case when 
one supposes only the fact of anisotropy. Using nonzero values of ct\ and a 2 one can 
also estimate the influence of these parameters on the borderline dimension d c . It is 
interesting to calculate the dependence of d c on the parameter e too. In Fig. |l], this 
dependence and the dependence on small values of ot\ and a 2 are presented. As one 
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can see from this figure d c increases when e — > and also the parameters a.\ and a 2 
give small corrections to d c . In Fig. ^, the dependence of d c on a x and «2 for e = 2 
is presented. 



4 Conclusion 

We have investigated the influence of the weak anisotropy on the fully developed 
turbulence using the quantum field RG double expansion method and introduced 
the modified minimal substraction scheme in which the space dimension dependence 
of the divergent parts of the Feynman diagrams is kept . We affirm that such 
a modified approach is correct when one needs to compute the d— dependence of 
the important quantities and is necessary for restoration of the stability of scaling 
regimes when one makes transition from dimension d = 2 to d = 3. We have derived 
analytical expressions for the fixed point in the limit of the weak anisotropy and 
found the equation which manages the stability of this point as a function of the 
parameters e, a% and a 2 , and allows one to calculate the borderline dimension d c . 
Below this dimension the fixed point is unstable. In the limit case of infinitesimally 
small anisotropy («i — > and a 2 — > 0) and in the energy pumping regime (e = 2) 
we have found the borderline dimension d c = 2.44. We have also investigated the 
e-dependence of d c for different values of the anisotropy parameters ai,a 2 and also 
the dependence of d c on the relatively small values of a.\ and a 2 for the physical 
value e = 2. 
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Appendix I 

The explicit form of the parameter A and functions a$, 6j, q, di and 6j (i = 1...4) for 



the divergent parts of diagrams 
1 



a 2 



2d(2 + d)(4 + d)(Q + d) 
x [-48 - 20d + 70d 2 + 30d 3 - 2ld 4 - 10rf 5 - d 6 

+ a 2 (24 + lGd - 22d 2 - 16d 3 - 2d 4 ) + a x (24 + 52d - Ad 2 - 50rf 3 - 20rf 4 - 2d 5 ) 
+ xi (-36 - 78d + 6d 2 + 75d 3 + 30rf 4 + 3d 5 ) 

+ X2(-36 - 24rf + 33d 2 + 24rf 3 + 3d 4 ) + Xs(-36 - 9d + 36d 2 + 9d 3 )] , 
1 

4rf(2 + d) (4 + d) (6 + d) 
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x [t*i(-96 - QAd + 88d 2 + 64d 3 + 8d 4 ) 

+ a 2 (-96 - QAd + I24d 2 + 82d 3 - 2Qd 4 - 18d 5 - 2d 6 ) 

+ xi(144 + 96d- 132d 2 -96d 3 - 12d 4 ) 

+ X2 (144 + 96d - 186d 2 - 123d 3 + 39d 4 + 27d 5 + 3d 6 ) 

+ xs(72 + Qd - 87d 2 - 9d 3 + 15d 4 + 3d 5 )} , 

a 3 = a 2 , 

6X3(1 -d 2 ) 

° 4 " (2 + c*)(6 + d)' 



d(2 + d)(4 + d)(6 + d) 
x [-48 - 20d + 70d 2 + 30rf 3 - 21d 4 - 10d 5 - d 6 + a 5 (12 + 3d - 12d 2 - 
+ (a 2 + a 4 ) (12 + 8d - 1 Id 2 - 8d 3 - d 4 ) 
+ («i + a 3 ) (12 + 26rf - 2d 2 - 25d 3 - 10rf 4 - d 5 ) 
+ xi (-36 - 78d + Qd 2 + 75d 3 + 30d 4 + 3d 5 ) 

+ x 2 (-36 - 24d + 33d 2 + 24d 3 + 3d 4 ) + % 3 (-36 - 9d + 36d 2 + 9d 3 )} , 



2d(2 + d)(A + d)(Q + d) 
x [(«! + a 3 ) (-48 - 32d + 44rf 2 + 32rf 3 + Ad 4 ) 
+ a 5 (-24 - 2d + 29rf 2 + 3d 3 - 5d 4 - d 5 ) 
+ (a 2 + «4) (-48 - 32d + 62d 2 + 41d 3 - 13d 4 - 9d 5 - d 6 ) 
+ xi(144 + 96rf- 132d 2 -96d 3 - 12rf 4 ) 
+ X2(144 + 96d - 186d 2 - 123d 3 + 39d 4 + 27d 5 + 3d 6 ) 
+ xs(72 + 6d - 87d 2 - 9d 3 + 15d 4 + 3d 5 )} , 
h = b 2 , 

4(rf 2 -l)(a 5 -3 X3 ) 
4 (2 + d)(Q + d) ' 

1 

Cl ~ 2d(2 + d)(A + d){Q + d) 

x [-48 - 20d + 70d 2 + 30d 3 - 21d 4 - 10d 5 - d 6 

+ a 5 (24 + Qd - 24d 2 - Qd 3 ) + a 4 (24 + lQd - 22d 2 - 16d 3 - 2d 4 ) 

+ a 3 (24 + 52d - 4d 2 - 50rf 3 - 20d 4 - 2d 5 ) 

+ xi (-36 - 78d + Qd 2 + 75d 3 + 30d 4 + 3d 5 ) 

+ X 2 (-36 - 24d + 33d 2 + 24d 3 + 3d 4 ) + Xs(-36 - 9d + 36d 2 + 9d 3 )} , 
1 

° 2 ~ 4d(2 + d)(4 + d)(Q + d) 

x [a 3 (-96 - 64d + 88d 2 + 64d 3 + 8d 4 ) 

+ a 5 (-48-4d + 58d 2 + 6d 3 -10d 4 -2d 5 ) 
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+ a 4 (-96 - 64d + I24d 2 + 82d 3 - 26d 4 - 18d 5 - 2d 6 ) 

+ xi(144 + 96rf- 132d 2 -96d 3 - 12d 4 ) 

+ X2 (144 + 96d - 186d 2 - 123d 3 + 39d 4 + 27d 5 + 3d 6 ) 

+ xs(72 + 6d - 87d 2 - 9d 3 + 15d A + 3d 5 )] , 

C3 = c 2 , 

(rf 2 -l)(4a 5 -6x 3 ) 

C4 ~ (2 + d)(6 + d) ' 



4d(2 + d)(4 + d)(6 + d) 
x [24d - 14d 2 - 33d 3 + 13d 4 + 9d 5 + d 6 + a 2 (12 - 4d - 13d 2 + 4d 3 + d A ) 
+ a!(-12-20rf + 3rf 2 + 19rf 3 + 9rf 4 + rf 5 ) 
+ xi (36 + 42d - 18d 2 - 40d 3 - 18d 4 - 2d 5 ) 
+ X2(-12 + 16d + 15d 2 - 16d 3 - 3d 4 ) + x 3 (6 + 9d - Qd 2 - 9d 3 )} , 



8d(2 + d)(4: + d)( y 6 + d) 
x (-48 + 16d + 52d 2 - 16d 3 - 4d 4 ) 
+ a 2 (48 + 80d - 60d 2 - 96d 3 + 10rf 4 + 16rf 5 + 2d 6 ) 
+ xi (48 - 64d - 60d 2 + 64d 3 + 12d 4 ) 
+ X2M8 - 104d + 62d 2 + 127d 3 - lid 4 - 23d 5 - 3d 6 ) 
+ xs(-2d + 7d 2 + 5d 3 - 7d 4 - 3d 5 )} , 



8d(2 + d)(4 + d)(6 + d) 
x [ai(48 + 56d - 40d 2 - 56d 3 - 8d 4 ) + a 2 (-48 - 56d + 40d 2 + 56d 3 + 8d 4 ) 
+ xi (-48 - 104d + 32d 2 + 104d 3 + 16d 4 ) 
+ X2(48 + 32d - 38d 2 - 25d 3 - 9d 4 - 7d 5 - d 6 ) 
+ xs(22d-d 2 -21d 3 + d 4 -d 5 )], 
X 3 (-10 + d + 10d 2 - d 3 ) 
2(2 + d)(6 + d) 
1 

4d(2 + d)(4 + d)(6 + d) 
x [24d - Ud 2 - 33d 3 + 13d 4 + 9d 5 + d 6 + 3da 5 (-l + d 2 ) 
+ a 4 (12-4d- 13d 2 + 4d 3 + d 4 ) 
+ a 3 (-12-20d + 3d 2 + 19d 3 + 9d 4 + d 5 ) 
+ xi (36 + 42d - 18d 2 - 40d 3 - 18d 4 - 2d 5 ) 

+ X2(-12 + 16d + 15d 2 - 16d 3 - 3d 4 ) + x 3 (6 + 9d - 6d 2 - 9d 3 )] , 
1 

8d(2 + d)(4 + d)(6 + d) 
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e 3 



e 4 
.4 



x [a 3 (-48 + 16d + 52c/ 2 - 16c/ 3 - 4c/ 4 ) + a 5 (-8c/ 2 - 2c/ 3 + 8c/ 4 + 2c/ 5 

+ a 4 (48 + 80c/ - 60c/ 2 - 96c/ 3 + 10c/ 4 + 16c/ 5 + 2c/ 6 ) 

+ xi(48 - 64c/-60c/ 2 + 64c/ 3 + 12c/ 4 ) 

+ x 2 (-48 - 104c/ + 62c/ 2 + 127c/ 3 - 11c/ 4 - 23c/ 5 - 3c/ 6 ) 

+ X3(-2c/ + 7c/ 2 + 5c/ 3 -7c/ 4 -3c/ 5 )], 
1 

8c/(2 + c/)(4 + c/)(6 + c/) 
x [24c/a 5 (-l + c/ 2 ) + a 3 (48 + 56c/ - 40c/ 2 - 56c/ 3 - 8c/ 4 ) 
+ a 4 (-48 - 56c/ + 40c/ 2 + 56c/ 3 + 8c/ 4 ) 
+ xi (-48 - 104c/ + 32c/ 2 + 104c/ 3 + 16c/ 4 ) 
+ X2(48 + 32c/ - 38c/ 2 - 25c/ 3 - 9c/ 4 - 7c/ 5 - c/ 6 ) 
+ X3(22c/-c/ 2 -21cZ 3 + cZ 4 -c/ 5 )], 

6a 5 (l - c/ 2 ) + X3(-10 + d + 10c/ 2 - c/ 3 ) 
2(2 + d)(6 + d) ' 



(2vr) d (c/ 2 - 1) ' 

where Sd is c/-dimensional sphere given by the following relation: 

27rf 



Sd 



r(i: 



Appendix II 

We present here the explicit analytical expressions for the fixed point in the weak 
anisotropy limit and also the equation which governs its stability. 
The basic form of the fixed point is 



9u 


= 010* + 9n*ai + gi2*a 2 


92* 


= 920* + g 2 l*a\ + 022*«2 


«3* 


= en«i + ei2« 2 , 


«4* 


= e 2 i«i + e 22 a 2 , 


Xi* 


= e 3 ioti + e 32 a 2 , 


X2* 


= e^oti + e 42 a 2 , 


«5* 


= o, 


X3* 


= o, 



where g w * and g 2 o* are defined in eq.([23|), and gn*,gn*, g 2 i*, g 2 2* and e^, i = 1,2, 3, 4, 
j — 1,2 are functions only of the dimension c/ and parameters e and 5 = (d — 2)/2. 
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They have the following form: 



giln _ 9l2n 921n _ 922n 

9n* — — ) 912* — — , g2i* — — , g22* — — , 
gud g\2d g22d g22d 

e lln _ &12n _ e 2ln _ &22n 
e ll — , Ci2 — , &21 — , ^22 — , 

ed e d e d e d 

e 31n _ C32n _ &41n _ &A2n 

e 3 i — , e 32 — , 641 — , e 42 — , 

g s ed g s ed g s ed g s &d 

where 

gun = 3(d 2 - l)gio*(d 6 (gio* + #2o*)(( 5e 3i - 3)510* - Se n g 2 o* + 5e 3 i5 2 o*) 

+ 3rf 5 (5( 10 * + 5 f 2o*)((— 2 + 3e 31 + 2e 41 )g w * - (2e n + e 21 - 3e 31 - 2e 4 i)5 2 o*) 

- 8(#i * + 5 , 2o*)((— 1 + 3e 3 i - e 4 i)^i * - (en - e 21 - 3e 31 + e 41 )g 2 o*) 
+ ^(-((5(10* + 5'20*)((-4 + 9e 3 i - 6e 4 i)5i * 

+ (-4eii + 3(e 2 i + 3e 3 i-2e 4 i))#20*)) 

+ 85((-5 + 10e 3 i + 3e 4 i)fi'io* - (5en + e 21 - 10e 31 - 3e 4 i)g 20 *)) 

+ 2d(-((g 10 * + ^20*) ((-1 + 6e 4 i)gio* - (en + 3e 2 i - 6e 4 i)y2o)) 

+ 165((-1 + 3e 34 - e 4 i)^i * - (e n - e 21 - 3e 31 + e 4 i)^ 20 *)) 

+ d 2 ((g w * + 5'20*)((— 15 + 34e 3i )5(io* + (-15en + 5e 2 i + 34e 31 )g 20 *) 

+ 165((-4 + 9e 31 + 2e 4 i)# 10 * 

+ (-4en + 9e 31 + 2e 4 i)5 2 o*)) + d 4 (8<5(-5i * + 2e 3 i5i * - en^o* + 2e 3 i5 2 o*) 

- (#io* + 02O*)((-1O + 15e 3i + 8e 4 i)5fi * 
+ (-10e n - 3e21 + 15e 31 + 8e 4 i)ff 2 o))) , 

g lld = 2d(4 + d)(-15d 6 (g 10:¥ + g 20 *) 2 + 6d 7 (g w * + g 20 *) 2 

+ 2(5(10* + 520*) (16e - 3(510* + 520*)) + 4rf 4 (e(5i * - 2520*) 

+ 6(510* + 520*) 2 + 35(25io* + 520*)) + d 5 (3(5io* + 520*) 2 + 125(25 10 * + 520*) 

- 4e(5 10 * + 4520*)) - 4rf 3 (6(5i * + 5 2 o*) 2 ~ 3e(5i * + 45 20 *) 

+ 5(8e + 9(25io* + 520*))) + d(15(5i * + 520*) 2 - 8e(5i * + 4520*) 

+ 5(-128e + 24(25io* + 520*))) - d 2 (A5(32e + 6510* + 35 20 *) 

+ 3((5i * + 520*) 2 + 4e(35i * + 2520*)))) , 

g2m = -((-1 + d 2 )(-3(-4-d + 4d 2 + d 3 )5io*(2(-2 + d 2 )5i * 

- (4 - 3d + d 2 )g 20 *) x 

x ((6e 3 i + d 2 (-l + 2e 3 i) - 2(1 + e 4 i) + 3d(-l + 2e 3 i + e 4 i))^io* 

- ((2 + 3d + d 2 )e u 

+ (-2 + d)e 21 - 6e 3 i - 6rfe 3 i - 2<i 2 e 3 i + 2e 4i - 3rfe 4 i)5 2 o*) 

+ d(4 + d)((-4 + 6e 3 i + d 2 {-2 + 3e 31 ) + 6e 4 i + d(-8 + 12e 31 + 3e 41 ))^ 2 , 

+ (2 + d 2 {l - 2en) - 4en - 4e 2 i - 6e 3 i + d(l - 8en - 2e 2 i + 6e 3 i - 3e 44 ) 

+ 18e 4 i)5i *52o* + ((2 + d + d 2 )e u + (-10 + d)e 21 - 3(4e 3 i + 2rfe 3 i + rf 2 e 3 i 
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- 4e 41 + 2de 4 i))(? 2 2 J(-8(2 + d)e + 3(-l + d)\l + d)(2g w * + g 20 *)))) , 

92id = 3(-l + d) 2 d(4 + 5d + d 2 )(-4-d + 4d 2 + d 3 )g 10 M-2 + d 2 )gw* 

- (4 - 3d + d 2 )g 20 *) + d{4 + d) 2 (-8(2 + d)e 

+ 3(-l + df(l + d)(2g 10 * + g 2 o*))(d 2 (85 + 3g 10 *) - 4(g 10 * + 520*) 

- 3d 3 (g Wif + 25(20*) + d 4 (g w * + 4g> 20 *) + d(16<5 + 35i * + 6520*)) , 

9\2n = 3(-l + C? 2 )fi'l0*(3rf 5 (5'l0* + 02O*)(( -1 + 3e 32 + 2e 42 )5i * 

- (2ei2 + e 22 - 3e 32 - 2e 42 )5 2 o*) - 8(510* + 5 , 2o*)((l + 3e 32 - 642)5*10* 

- (ei2 - e 22 - 3e 32 + e 42 )5 2 o*) 

+ d 6 (g 10 * + 520*) (-3612520* + 5e 32 (5io* + 520*)) 

+ rf 3 (-((5i * + 520*)(3(1 + 3e 32 - 2e 42 )5i * 

+ (-4ei2 + 3(e 22 + 3e 32 - 2e 42 ))5 2 o*)) 

+ 8S((-1 + 10e 32 + 3e 42 )5i * - (5e 12 + e 22 - 10e 32 - 3e 42 )5 2 o*)) 

+ 2d(— ((^10* + #20*) ((-3 + 6e 4 2)5i * - (e 12 + 3e 22 - 6e 42 )5 2 o*)) 

+ 165((1 + 3e 32 - e 4 2)5i * - (e 12 - e 22 - 3e 32 + 642)520*)) 

+ d 4 (-((5i * + 520*) ((-3 + 15e 32 + 8e 42 )5i * 

+ (-10ei2 - 3e 22 + 15e 32 + 8e 42 )520*)) 

+ 85(-(ei2520*) + 2e 32 (5i * + 5-20*))) + d 2 ((g w * + 520*) ((5 + 34e 32 )5i * 

+ (-15ei2 + 5e 22 + 34e 32 )5 2 o*) + 165(9e 32 (5i * + 520*) 

+ 2 (-2ei252o* + e 4 2 (910* + #20* ))))), 

9l2d — 9lld j 

922n = -((-l + rf 2 )(-3(-4-rf + 4rf 2 + rf 3 )5io*(2(-2 + rf 2 )5io* 

- (4 - 3d + d 2 )52o*)((2 + 6e 32 + 2d 2 e 32 - 2e 42 + d(-l + 6e 32 + 3e 42 ))5io* 

- ((2 + 3d + d 2 )e 12 + (-2 + d)e 22 

- 6e 32 - Qde 32 - 2d 2 e 32 + 2e 4 2 - 3de A2 )g 20 *) 

+ d(4 + d)((-4 + 6e 32 + 3d 2 e 32 + 6e 42 + d(-2 + 12e 32 + 3e 42 ))5i 2 o* 

- (10 + 4ei2 + 2rf 2 ei2 + 4e 22 + 6e 32 - 18e 42 
+ d(-l + 8ei2 + 2e 2 2 - 6e 32 + 3e 4 2))5io*fi'2o* 

+ ((2 + d + d 2 )e 12 + (-10 + d)e 22 - 3(4e 32 + 2cfe 32 + d 2 e 32 

- 4e 42 + 2rfe 42 ))52 2 o*)(-8(2 + d)e + 3(-l + d) 2 {\ + d){2g w * + 520*)))) , 

#22^ = 92\d , 

eiin = (g q - gsP2)(g P g s (m 4 n 2 -m 3 n 3 )p 1 

+ +gio*g ((iTi 4 n 1 + min 3 )p 4 - (m 3 rai + m 1 n 2 )p 5 )) , 

e d = g 3 s (m 4 n 2 - m 3 n 3 )pip 2 + g20*g o g q (-(rn A n 1 + min 3 )p 4 + (m 3 ni + min 2 )p b ) 

+ g20*g o gs((m 1 n 3 p 2 + m 4 rai(-pi + p 2 ))p 4 + m 3 nipip 5 - m 3 n x p 2 p^ 

- m x n 2 p 2 p h + m 2 p x {n 3 p A - n 2 p 5 )) + gk{g q g s {-{m A n 2 ) + m 3 n 3 ) 
+ g2Q*g {m A nip A - m 2 n 3 p A - m 3 nip 5 + m 2 n 2 p 5 )) , 
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ei2n = (g q ~ gsP2)(g P gs(m 4 n 2 -m 3 n 3 )p 2 

+ +gio*g (-m4n 1 p 4 + m 2 n 3 p 4 + m 3 n x p^ - m 2 n 2 p 5 )) , 

e 2 in = (gk - gsPi)(g P g s (m 4 n 2 - m 3 n 3 )pi 

+ +gio*g {m A n l p A + m x n 3 p A - m 3 riip^ - min 2 p 5 )) , 

e 2 2n = {gk - g s Pi)(g P g s (m 4 n 2 - m 3 n 3 )p 2 

+ gio*g (-m4n 1 p 4 + m 2 n 3 p 4 + m 3 n x p h - m 2 n 2 p 5 )) , 

e3in = g2Q*g P g s Pi{-gkm A n 1 + g q m A n x + g q m x n 3 + g k m 2 n 3 + g s m A nipi 

- g s m 2 n 3 p x - g s m A n x p 2 - g s m x n 3 p 2 ) + gw*(g s Pi(gs( m 4 n i + ™i™3)p 2 

- g2Q*g (m 1 + m 2 )nip b ) + ft(-(3,9 s (m 4 ni + TOin 3 )) 
+ g20*g o (m 1 + m 2 )nip 5 )) , 

e32n = g2o*g P g s P2(-gkm A n 1 + g q m A n x + g q m x n 3 + g k m 2 n 3 + g s m A n x p\ 

- g s rn 2 n 3 px - g s m A nYP2 - g s min 3 p 2 ) + gio*(gkg q g s (rn 4 ni - m 2 n 3 ) 
+ m 4 ri! + m 2 n 3 )p x p 2 - g 2 o*g g q (m 1 + m 2 )n x p^ 

+ g20*gog s (m 1 + m 2 )nip 2 p 5 ) , 

e 4 in = g2o*g P g s Pi(-gkm 3 n 1 + g q m 3 n x + g q m x n 2 + g k m 2 n 2 + g s m 3 nipi 

- g s m 2 n 2 p x - g s m 3 n x p 2 - g s m 1 n 2 p 2 ) + g w *{g s Pi{g 2 s {m 3 ni + m x n 2 )p 2 

- g2Q*g (mi + m 2 )nip A ) + g k (-(g q g s (m 3 n 1 + min 2 )) 
+ g2o*g (m 1 + m 2 )nip 4 )) , 

e 42 n = g2o*g P g s P2(-gkm 3 n 1 + g q m 3 n x + g q min 2 + g k m 2 n 2 + g s m 3 nipi 

- g s Tn2n2Pi - g s m 3 n x p 2 - g s m x n 2 p 2 ) + #10* (^^(m^i - m 2 n 2 ) 
+ gl(-m 3 ni + m 2 n 2 )pip 2 - g20*g o g q (m 1 + m 2 )nip A 

+ 920*Ms( m i + ra 2 )rciP 2 P4) , 

where 





= 24 + 16d - 22d 2 - 16d 3 - 2d 4 , 


mi 


= 48 - 16d - 52d 2 + 16d 3 + 4d 4 , 


m 2 


= -48 - 80d + 60d 2 + 96d 3 - 10d 4 - 16d 5 - 2d 6 , 


m 3 


-48 + 112d + 32d 2 - 130d 3 + 14d 4 + 18d 5 + 2d 6 , 


m 4 


= 48 + 104d - 62d 2 - 127d 3 + lid 4 + 23d 5 + 3d 6 , 


ni 


= 48 + 56d - 40d 2 - 56d 3 - 8d 4 , 


™2 


= 48 + 104d - 32d 2 - 104d 3 - 16d 4 , 


n 3 


= -48 + 16d + 10d 2 - 41d 3 + 35d 4 + 25d 5 + 3d 6 , 


Ol 


= 26 - 7d - 27d 2 + 7d 3 + d 4 , 


o 2 


= -12 + 12d 2 , 


Pi 


= -96 - 64d + 88d 2 + 64d 3 + 8d 4 , 




= -96 - 64d + 124d 2 + 82d 3 - 26d 4 - 18d 5 - 2d 6 , 
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P3 


= 96 + 40d - UOd 2 - 60d 3 + 42rf 4 + 20d 5 + 2d 6 , 


Pi 


= 144 + 96d - 132c/ 2 - 96d 3 - l2d A , 


P5 


= 144 + 96d - 186cZ 2 - 123d 3 + 39d 4 + 27d 5 + 3d 6 , 


P6 


= -24d - 52d 2 + Ad 3 + 50d 4 + 20d 5 + 2d 6 , 


Pi 


= 96 + 16d - 192d 2 - 56d 3 + 92d A + 40d 5 + Ad 6 , 


Qi 


= 96 + 16d - 156d 2 - 38d 3 + 58d 4 + 22d 5 + 2d 6 , 


T\ 


= 24 - Ad - 3Qd 2 + 2d 3 + 12d 4 + 2d 5 , 


r 2 


= 12 + 2d- 18d 2 - 3d 3 + 6d 4 + d 5 , 


rs 


= 12 - 6d - 18d 2 + 5d 3 + 6d 4 + d 5 , 


9s 


— 9io* + 920*, 


9 P 


— 9w* + 5 'io*/ 920*, 


9k 


= {glo*Pz) 1 920* + 920*P6 + #10*^7, 


9 Q 


= (-(dgloJi) + 9io*(9io*P3 + 92o*qi))/92o*, 


9o 


= 9s/ 920* ■ 



Stability of the fixed point is determined by the (2x2) block of the stability matrix 
which corresponds to /3-functions of a 5 and xz- The eigenvalue which responds for 
instability has the form: 

A = A + Ai«i + X 2 a 2 , 

dg 2 o*(gio* + 920*)oi ~ Vh + 9io*920*ri + g 2 *r 2 + gl^Jz 
8d(l2 + 8d 2 + d 2 )g 20 * 

A In 

a7' 

= dgj 0m (g lu + ^2i*)oiaAi + 920* (~*2 + gu*Vh. (g2o*n + 2g w *r 2 )) 

+ g2i*Vh (Vti - g\o* r 2 + gloJ?) , 

= 8d(12 + 8d + d 2 )g 2 20 ,Vh, 

X 2n 

a7' 

= dgl Qir (gi2* + g22*)oiVh + #20* (-*3 + gi2*Vh (g2o*n + 2#i *r 2 )) 
+ g22*Vh (Vh - g 2 1Q j2 + gloJz) , 

*i = d 2 gl Qif {gi^ + #2o*) 2 (°i - 4o 2 ) - 2dg 20 *o 1 (g 10 * + 5-20*) x 

X (#10*#20*?"l + 9lo* r 2 + ^20* r 3) + (#10*#20*?"l + g\o* r 2 + gloJzf , 



where 

Ao 
Ai 

Aln 

A d 
A 2 

A2n 

with 
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t 2 = 2{d 2 g2 *(gio* + 52o*)(5n*52o* + (#10* + 2# 2 o*)#2i*)(oi - ^o 2 2 ) 

+ (510*520*^1 + 5io* r 2 + 92o* r 3)(9u*92o*ri + 5io*52i*?"i + 2^10*511*^2 

+ 25-20*521*^3) - ^ol(5?o*52i*^2 + 5? *52o* (35n*^2 + 2521* ij\ + r 2 )) 

+ 52o*( 4 52i*r 3 + 011* (ri + r 3 )) + 510*520* ( 2 5n*(n + r 2 ) + 3# 2 i*(n + r 3 )))) , 

t 3 = 2(d 2 5-20*(5l0* + 520*) (512*520* + (510* + 2#20*)522*)(Oi - 402) 

+ (510*520*^1 + 5?0* r 2 + 5 20 * r 3) (512*520*^1 + 510*522*^1 + 2^10*512*^2 

+ 25-20*522*?" 3 ) - doi(g\ Qi ,g-22*r2 + 5io*520*(35i2*r 2 + 2g 22 ^n + r 2 )) 

+ 52o*( 4 522*?" 3 + 5i2*(n + r 3 )) + 510*520* ( 2 5i2*(?"i + r 2 ) + 3g22*( r i + r s)))) ■ 

borderline dimension d c is defined as a solution of the equation 

X(d c , e, ai, a 2 ) = . 
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